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1.  TECHNICAL  PROBLEM 


We  are  investigating  the  consequences  of  a  new  theory  of  gravita¬ 
tion  by  employing  the  symboMc  manipulation  system  MACSYMA  developed 
by  the  Mathematical  Laboratory,  Project  MAC,  MIT.  The  new  theory  was 
formulated  by  Dr.  Huseyin  Yilmaz  (See  Section  2)  over  the  last  twenty 
years  and  is  believed  to  represent  the  or1**  viable  alternative  to 
General  Relativity,  Einstein's  gravitational  theory. 

There  are  three  main  objectives  in  the  contract  which  are  currently 
being  investigated  using  MACSYMA.  The  theoretical  question  we  shall 
answer  concerns  the  verification  that  a  conjectured  relation  satisfies 
fMpId  pnnaH  nf  rho  fhnnry  Cr.  pC  cor-nnH  o^pr.  J.TnH_p.n 

the  usual  field  theory  in  which  field  equations  are  written  down  and 
subsequently  solved,  the  new  theory  presents  a  conjectured  solution 
(Reference  5,  Section  2,  or  Appendix  1)  of  the  field  equations  to  be 
verified  in  detail.  A  related  theoretical  question  which  we  hope  to 
answer  is  that  the  equations  of  motion  of  objects  follow  from  the  field 
equations  and  conservation  laws.  This  feature  has  never  been  well  im¬ 
plemented  in  the  Finsteinian  theory. 

The  new  theory  differs  from  General  Relativity  in  many  fundamento1 
ways.  We  have  shown,  for  example,  that  certain  types  of  gravitational 
radiation  believed  to  exist  in  General  Relativity  in  fact  are  absent  in 


that  theory  (Appendix  2).  Also,  it  seems  clear  in  Einstein's  theory 
that  gravita'  »  •  *•  waves  (if  they  exist)  cannot  carry  energy.  On  the 
other  hand  we  have  shown  that  gravitational  radiation  does  exist  in 
the  New  Theory  and  these  waves  do  carry  energy.  In  addition  we  are 
considering  the  New  Theory  as  a  unified  field  theory  in  which  gravita¬ 
tion  and/or  electromagnetism  arise  from  the  field  equations.  It  is 
hoped  that  these  studies  will  answer  the  question  of  whether  informa¬ 
tion  may  be  transmitted  via  gravitational  waves  and  whether  gravitation 
modifies  electromagnetic  rrdiation  to  any  practically  significant  extent. 


2.  LITERATURE  REVIEW 

The  basic  literature  concerning  the  new  theory  consists  essentially 
ot  the  lollowing  11  publications  by  Dr.  Huseyin  Yilmaz. 

1.  H.  Yilmaz,  Phys.  Rev.,  Ill,  p.  1417  (1958). 

2.  H.  Yilmaz,  "Evidence  for  Gravitational  Theories",  ed.  C.  Moller, 
Academic  Prars,  New  York  (1962). 

3.  H.  Yilmaz,  "Introduction  to  the  Theory  of  Relativity",  Blaisdell 
Publishing  Company,  New  York  (1965). 

4.  H.  Yilmaz,  Proceedings  of  the  Rochester  Meeting  of  the  Division 
of  Particles  and  Helds  of  the  American  Physiral  Society,  ed.  A.  C. 
Melissinos  and  P.  F.  Slattery,  America!  Institute  of  Physics,  p.  259,  1971. 

5.  II.  Yilmaz,  "New  Theory  of  Gravitation",  Physical  Review  Letters, 


vol.  27,  p.  1399,  November  1971. 


6.  H.  Yilmaz ,  Nuovo  Cimento,  10B,  p.  79,  1972. 


7.  H.  Yilmaz,  "On  the  New  Theory  of  Gravitation",  Lettere  A1 
Nuovo  Cimento,  vol.  5,  no.  4,  September  1972. 

8.  H.  Yilmaz,  "Correspondence  Paradox  in  General  Relativity", 
Lettere  A1  Nuovo  Cimento,  vo  1.  7,  no.  9,  p.  337,  June  1973. 

9.  H.  Yilmaz,  "Postulational  Derivation  cf  the  New  Theory  of 
Gravitation",  Lettere  A1  Nuovo  Cimento,  vol.  6,  no.  5,  p.  181,  Feb. 
1973. 

10.  H.  Yilmaz,  Physics  Today,  p.  15,  March  1973. 

11.  H.  Yilmaz,  Annals  of  Physics,  vol.  81,  pp.  179-200,  1973. 


3.  TECHNICAL  RESULTS  TO  NOVEMBER  30,  1973 


In  the  following,  we  present  an  outline  of  our  progress  to  date  in¬ 
cluding  new  results  we  have  obtained  as  well  as  the  new  capabilities  of 
the  MACSYMA  system  we  have  introduced.  We  have  shown  that  plane  gravi¬ 
tational  waves  in  vacuum  (contrary  to  current  belief)  do  not  exist  in 
Einstein's  theory  (Appendix  2).  Previous  calculations  proved  the  exis¬ 
tence  of  the  waves  in  Einstein's  theory  in  fr'-st  order  cnly.  Claims 
for  exact  plane  wave  sclutions  were  advanced  but  never  rigorously  demon¬ 
strated.  Our  surprising  result  that  such  waves  do  not  exist  in  Linstein' 
theory  follows  by  computing  the  Einstein  tensor  to  second  order  and 
beyond  using  the  MACSYMA  system. 


In  a  si.’ll  r  -•!/  it  is  shown  that  such  waves  do  indeed  exist  in 
the  New  Theory  .Appendix  2,  and  Reference  11). 

Using  the  MAOSYMa  system  we  have  further  found  the  general  static 
solutions  of  the  Einstein  vacuum  equations  and  Einstein’s  equations  with 
a  cosmical  term  in  spherical  coordinates  (Appendix  3).  We  have  also 
found  the  general  solution  of  Einstein's  vacuum  equations  in  Cartesian-like 
coordinates  (Appendix  4).  It  is  doubtful  that  such  solutions  could  be 
found  by  hand  calculations.  The  enormous  value  of  the  MACSYMA  system 
is  here  evident  in  helping  us  to  formulate  conjectures  and  Making  it 
possible  to  actually  test  them.  The  few  previously  known  solutions  are 
shown  to  be  special  cases  of  the  general  solution  for  appropriate  choices 
of  a  free  function  which  appears  in  the  general  solution.  In  addition, 
while  investigating  the  Eincteinian  theory  we  have  discovered  solutions 
which  seem  to  violate  the  very  foundations  of  the  theory  itself 
(Appendix  5). 

We  have  implemented  a  first  step  in  the  MACSYMA  system  to  allow 
indicial  symbolic  manipulation  (Appendix  6).  Using  these  new  programs 
we  have  verified  the  first  order  functional  expansion  described  above. 

This  is  a  major  development  from  the  standpoint  of  symbolic  nipulation 
(to  our  knowledge  no  other  system  has  this  facility).  We  shall  gradually 
expand  and  improve  this  program. 

We  have  made  some  progress  in  the  understanding  of  matrices  gene¬ 
rated  by  expansions  of  the  form  e*  =  1  +  X  +  ?  . 

<_  i  T  A  ^  2!  +  *  •*  whGre  *  is  a  square 


^trlx  (Appendix  7).  Such  considerations  are  necessary  for  the  study 

of  unified  field  theory  as  well  as  the  new  theory  In  its  present  form. 

Expansions  of  this  type  are  easily  generated  by  the  MACSYMA  system  and 

has  allowed  us  to  obtain  many  closed  forms  of  the  series  for  particular 
choices  of  X. 

We  have  begun  to  Investigate  stress-energy  tensors  generated  by 

metrics  with  skew-symmetric  parts.  This  is  part  of  the  effort  directed 
toward  a  unified  field  theory. 

the  follow!-.*  list  of  capabilities  of  the  MACSYMA  system  is  directly 
attributable  to  the  combined  effort  of  Perception  Technology  Corporation 
and  our  program  consultants  at  Project  MAC.  At  the  start  of  our  contract 
the  MACSYMA  system  was  only  capable  of  computing,  for  a  given  metric, 
the  Inverse  metric,  Christoffel  symbols  and  kicci  tensors.  The  latter 
of  these  had  a  bug  which  was  discovered  and  corrected  by  Perception 

Technology  Corporation.  The  list  of  capabilities  as  of  November  30, 

1973  include. 

1)  For  any  metric  tensor  the  Ricci  tensors  and  Einstein  tensors 

(covariant  or  mixed)  can  be  computed  In  CRE  (canonically  rational  ex- 
pression)  form. 

2)  For  any  metric  tensor,  the  Inverse  metric,  Christoffel  symbols 
ot  born  kinds,  and  mixed  Einstein  tensors  can  be  computed  in  CRE  or  in 
«  truncated  power  series  to  arbitrary  order  in  an  expansion  parameter. 

We  are  limited  only  in  the  memory  of  the  system. 


3)  For  any  metric  tensor  the  compoi  ^r.-  <>'  the  full  Riemann  tensor 
can  be  computed  in  CRE  or  in  a  truncated  po- .  .  varies. 

A)  For  any  metric  tensor  the  equations  oL  im-;.1  nn  can  computed  in 
CRE  or  in  a  truncated  power  series. 

5)  For  any  metric  tensor  the  d ' Alemhertian  and  the  tV  of  the  new 

V 

theory,  both  applied  to  a  scalar  field,  can  be  computed. 

6)  For  any  metric  tensor  the  components  of  any  second  rank  tensor 
in  a  transformed  coordinate  system  can  be  computed. 

7)  For  a  function  <*>  (z-t)  one  has  the  differential  identities 

-  “  tfv.  $  =  -  <f>  ,  4>  =  +  <J>  which  are  the  wave  conditions  of 

2  L  ZL  LL  ZZ  tt 

certain  specially  chosen  metrics .  For  components  of  the  metric  tensor 
containing  functions  cf  this  kind  the  wave  conditions  may  now  be  auto¬ 
matically  implemented  with  MACSYMA  thus  trpmendous]v  shortening  compu¬ 
tations  concerned  with  gravitational  radiation. 

8)  The  indicial  manipulation  program  is  now  able  to  compute  all 
Ricci  and  Einstein  tensor  components  in  first  order.  Soon  the  capa¬ 
bility  will  extend  to  second  order. 


7 


Volumi  27,  Ni'muh)  20  PHYSICAL  R  I  V  1 1  \\r  LETTERS 


1 5  Novfmbih  1971 


APPENDIX  1 

New  Theory  of  Gravitation 

Hiiseyin  Yilmaz 

Perception  Technology  Corporation,  Winchester,  Massachusetts  01S90 
(Received  2  June  J  071 ) 


A  W-nIJ  V  I.OTnnt?-hM"tHonf  f'M^OfLrn^Ap  tknewv  esf  ,  ,vo-.  Al-  I  ,  f,  ,  ,  , 

_  |-W  ’  J  •  FI'V  *  v/C  tu  1IV1U  lO 

a  massless,  spin-2  field  vlt  of  PnuJi-Ficr?.  tvjx*  is  presented.  In  this  theory  is  no 
longer  the  gravitational  potential  itself  but  reduces  to  a  functional  gm.(Va1')  of  y>  Sev¬ 
eral  rigorous  and  a  general  iteration  solution  of  this  type  arc  exhibited.  The  static  cen¬ 
tral  bcxly  problem  reduces  exactly  to  author’s  1953  theory  in  the  special  case  <p -V>0° 

-  <p  and  a  (gf,  b  <p  -  o  so  that  as  in  that  theory  the  three  crucial  tests  are  satisfied. 


About  ten  years  ago  there  existed  at  least  sev¬ 
en,  more  or  less  plausible,  theories  of  gravita¬ 
tion.  Three  of  these  were  the  flat-space  theories 
of  Poincare',  Whitehead,  and  Birkhoff.  In  the 
curved-space  category  were  the  theories  of  Ein¬ 
stein,  Jordan,  Brans-Dicke,  and  the  author,  Ein¬ 
stein’s  theory  being  of  course  the  most  prominent 
among  them.  In  1961  with  the  measurement  of 
gravitational  ted  shift  in  earthbound  laboratories 
it  became  clear  that  the  flat-space  theories  have 
to  be  abandoned.1  This  situation  led,  since  1961, 
to  a  vigorous  development  of  Jordan  and  Dicke 
type  theories  (Brans-Dicke  theory  in  fact  belongs 
to  a  class  of  Jordan  theories)  and  until  recently^ 
Brans-Dicke  theory  was  considered  quite  serious¬ 
ly  in  many  quarters. 

In  the  fall  of  1970,  however,  the  experiment  on 
the  time  delay  of  radar  sign. .Is  reflected  from 
Venus  and  Mercury  reached  accuracies  high 
enough  to  permit  at  least  a  preliminary  judge¬ 
ment  against  the  Jordan  and  Brans-Dicke  type  of 


theory.2  If  the  present  trend  of  results  continues, 
the  1916  theory  of  Einstein  and  the  largely  unde¬ 
veloped  1958  theory  of  the  presem  author  will 
probably  be  the  only  two  generally  known  and  ex¬ 
perimentally  viable  theories  of  gravitation.3  On 
the  other  hand  the  1958  theory  of  the  author,  as 
published,  has  the  theoretically  objectionable  fea¬ 
ture  that  it  does  not  possess  time-dependent  dy¬ 
namic  solutions.  This  of  course  casts  rerious 
doubts  on  the  validity  of  this  theory  leaving  Ein¬ 
stein’s  1916  theory  practically  unrivaled.  We  are 
therefore  motivated  to  present  a  natural  extension 
of  our  1958  theory  which  seems  to  open  a  new  and 
interesting  avenue  of  inquiry  as  well  as  avoiding 
the  above  mentioned  objection. 

Formulation  of  the  theory, — One  of  the  main 
differences  between  Einstein’s  theory  and  the 
present  theory  of  gravitation  is  that  in  Einstein’s 
theory  the  components  of  the  metric  tensor  g 
are  considered  as  the  components  of  the  gravita¬ 
tional  potential  and  are  therefore  functions  of 
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space-time  variables  as /;  (x),  whereas  in  the 
nev.  theory  the  gravitational  potentials  form  a 
different  field  c  J  and  /, ,  are  functions  or  func  ¬ 
tionals  of  these  as  £  ).  Thus  in  a  sense  the 

new  approach  considers  gravitation  as  a  local 
field  theory  in  a  curved  geometry  i  whose  cur¬ 
vature  arises  from  the  existence  of  the  fields 
themselves. 

With  this  conception  in  mine  it  is  possible  to 
formulate  the  theory  of  gravitation  directly  as  a 
functional  relation  between  the  fields  c/  and  the 
metric  .  We  shall  state  this  relation  in  the 
form  of  a  functional  differential  equation, 

<4’„,  =2 U\.dc  -£vatt<rl°-xlo*/c>utr),  (l) 

where  c  is  the  trace  of  c  t,' .  Equa  ion  (1)  is  here 
intended  to  represent  a  i  elation  between  the  com¬ 
ponents  of  ar.d s'  in  c  ertain  preferred  coor¬ 
dinate  system  and  therefore  it  is  not  a  priori  co¬ 
ordinate  invariant.  The  preferred  coordinates 
thus  introduced  will  1."  seen  to  have  an  interpre¬ 
tation  as  gauge  fields  satisfying  a  sourceless 
wave  equation.  In  accordance  with  our  local  re¬ 
quirements  above,  the  ;; , ,  may  be  asscined  to  re¬ 
duce  to  the  special  relativistic  form  r,  ,  =  (l,  -  i, 
-1,-1)  when  (U  /  -C.  In  other  words  in  the  ab¬ 
sence  of  the  fields  or  their  variations  the  metric 

it;  p]CC|iinArJ  ♦  n  rrm  m-c*.  J-x* ,  ■  ~  T  r  , 

^  --  ..  — w..x.  id  4ii  milui 

can  be  clone  without  loss  of  generality.’  The 
fields  c,/  themselves  will  then  arise  via  the  dis¬ 
tributions  of  matter  according  to  the  locally  spe¬ 
cial  relativistic  equation.. 

<V  =  /(<»«„»  ‘  )’dV'/r',  (2) 

where  the  prime  implies  the  usual  retardation 
consideration  on  the  integral.  We  shall  therefore 
have  the  propagation  equations 

□V,/'  *450"^',  (3) 

where  on  u'  is  the  usual  “matter”  tensor  and  r  2 
is  the  general  d’Alemberti  in  L2  =(v  -a)'1  3, (/-*,' 
X-C'"  3 ,)  which  under  (7)  below  turns  out  to  be  C  2 
i  =3‘01.  The  symbol  3(.  denotes  ordinary 
derivative  in  the  coordin:  tes  in  which  (1)  holds. 

'Vc  may  now  point  out  t  lat  in  such  a  gravitation¬ 
al  field  the  energy-momentum  “vector”  of  matter 
alone  cannot  bo  a  conserved  quantity  and  in  fact 
it  cannot  be  a  true  vector.  We  must  ado  a  “field” 
stress-energy  term  /  loanin'  in  order  to  ob¬ 
tain  a  conserved  energy-momentum  vector.  We 
write 

r„"  =ouy  Ml/4r)V,  (4) 

where  /,/  is  net  possible  to  determine  uniquely 


from  the  usual  special  relativistic  field  theo;y 
arguments  because  the  divergence  of  an  arbiir 
vector  can  be  added  to  the  Lagrangian  and  this 
wil>  yield  a  different  although  equivalent  expres¬ 
sion.  We  shall,  however,  require  that,  as  com¬ 
puted  directly  through  (1),  the  geometric,  diver¬ 
gence-free  tensor  R. '  -  J  R  will  be  equal  to 
8n7  „* .  In  other  words  we  require  the  field  equa¬ 
tions  to  be  of  the  form 

V  - '■Vk=2(c\v+V')=8*V  .  (5) 

In  this  way  we  are  able  to  guarantee  the  eonserva- 
tion  of  the  energy-momentum  vector  identically 
and  at  the  same  time  determine  the  form  of  t ' 
essentially  up  to  an  additive  cosmological  con¬ 
stant  A.  We  note  again  that  under  the  combined 
effect  of  (1)  and  (7)  the  covariant  derivative  of  the 
left-hand  side  of  (5)  is  expressed  in  terms  of  the 
ordinary  derivatives  to  yield  the  middle  expres¬ 
sion  so  that  the  conservation  can  hold  in  (5)  with 
ordinary  rather  than  covariant  derivatives  be¬ 
cause  of  (1)  and  (7).  In  fact  under  (7)  below  the 
empty -space  value  of  R  (o=0)  gives  the  Lagran- 
gim, 

Li  *  f  =  “  5  V  o*9  >  «  °  +  i  3  Va  v  t’  (6) 


yielding,  indeed,  a  familiar  special  relativistic 
V  corresponding  to  the  fields  c/,  <p. 

It  is  to  be  noted  carefully  that  in  order  for  the 
second-order  derivatives  of  </>„'  to  combine  prop¬ 
erly  into  !<,->  •  the  -lcld  must  satisfy  a  subsid¬ 
iary  condition 

a\V  =0,  (7) 

which,  contrary  to  a  first  impression,  does  not 
imply  a  separate  conservation  law  for  the  matter 
part  alone.  It  rather  implies  that  locally  the 
grav/ational  field  will  propagate  with  the  velocity 
of  light  c  in  all  directions.  Since  the  absolute 
values  of  C3,/  arc  unobservable  and  since  under 
(1)  and  (7)  C  2  no  longer  contains  quadratic  first- 
order  derivatives  of  </>/,  the  implied  local  coor¬ 
dinates  are  interpretable  as  sourceless  gauge 
fields  quite  analogously  to  special  relativity 
where  v,  y,  z,  and  t  satisfy  a  sourceless  wave 
equation.  It  is,  however,  essential  to  keep  (7)  as 
a  separate  condition  so  that  it  can  be  dropped  to 
allow  general  coordinate  transformations. 

finally,  we  adopt  the  geodesic  equations  of  mo¬ 
tion 


(P 


=  0 


(8) 
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as  our  equations  governing  the  motions  of  test 
particles  and  light  r  ...  ible  to  do  this 

because  any  line  elenu  •  i  :s  solution  of  (1) 
and  (7)  is,  by  construct ,  •  -o  a  solution  of  (5). 
Since  under  certain  gorcr-.I  conditions  the  geode¬ 
sic  equations  above  are  know:  o  L*-  a  consequence 
of  field  equations  (5)  we  are  jus'ifnd,  as  in  Ein¬ 
stein’s  theory,  in  postulating  (8;  as  our  equations 
of  motion.  (See  also  our  concluding  remarks  at 
the  end  of  the  paper.)  The  theory  thus  introduced 
goes  over  into  our  1958  form  in  the  limit  of  static 
fields  0h’  -  <pn°  -  tf  ,  d V  =  0. 

Comj'crisoH  with  Einstein's  theory. — It  is  now 
possible  to  start  comparing  the  above  theory  with 
Einstein’s  theory  to  see  that  the  two  theories  are 
very  different  conceptually  and  yet  they  give  the 
same  answers  for  the  three  classical  tests  of  gen¬ 
eral  relativity.  For  the  time  delay  of  radar  sig¬ 
nals  echoed  from  Venus  and  also  for  the  experi¬ 
ment  on  the  isotropy  of  inertia  there  seem  to  be 
certain  differences  in  prediction  or  interpreta¬ 
tion.  The  theoretical  differences  arc  mainly  re¬ 
lated  to  tiie  fact  that  in  Einstein’s  theory  one  as¬ 
sumes  the  right-hand  side  of  (5)  to  be  the  “mat¬ 
ter”  tensor  alone,  hence  t  u‘  does  not  enter.  Con¬ 
sequently  in  empty  space  where  cr=0  the  matter 
tensor  vanishes  as  well  and  one  has  the  Einstein- 
tan  law 

=°-  o) 

By  solving  these  equations  for  the  static  gravi¬ 
tational  field  of  the  sun  cue  obtains  the  usual 
Schwarzschild  line  elements.  In  the  present  the¬ 
ory  this  static  case  corresponds  via  (3)  to  y0n=  o 
=  M/r,  and  putting  this  into  (1)  we  immediately 


h’n  SV„,  +2if)((1,<P  -  +  4/j?;(J|  <p  - 

-  4  / 1 *  V  -  <P  X  “ -  lha  <?  ,  °]  dtp,  x 


have  the  line  element 

ds2  =e  '2^dt2  - r2Vv2  +d\2  +dz2),  (10) 

The  reader  could  check  easily  that  upon  comput¬ 
ing  the  tensor  li p’  —  {6^'R  with  this  line  element 
one  obtains  for  t ^  exactly  the  expression  (6)  with 
~  v0°  =  tp.  He  could  also  •  apeat  this  process 
for  other  solutions  given  below.  We  can  see  that 
the  two  theories  are  in  fact  \  ex j  different.  For 
example,  in  our  theory  we  now  have,  instead  of 
R/  =0,  the  new  empty  space  law 

V=2  (V'-ifi/O  (11) 

exhibiting  a  local  field  theory  form  explicitly. 

We  may  point  out  immediately  that  we  do  not 
have  to  solve  the  highly  nonlinear  equations  (11) 
directly  in  order  to  obtain  a  line  element.  In¬ 
stead  we  can  simply  deduce  it  from  our  functional 
equations  (1)  as  we  have  just  done  to  obtain  the 
static  line  element  (10).  As  other  examples  of 
such  exact  solutions  let  all  components  of  y>p" 
other  than  <f>2(t,  z)  vanish.  Then  the  solution 
is  A'oo  =  -£33  *  1.  ~S 22  =  ~  cosh4^2,  gn  =£21 

-sinh4</’,2.  If  only  <P03\v,y)  exists  then^,,  =gn 
=  “1.  ^rooi:--er33  =  COs4^03.  =S>n4‘iP03.  etC. 

The  latter  solution  with  a  slight  extension  is  found 
to  provide  an  explanation  (through  an  argument 
similar  to  that  given  by  Lense  and  Thirring  in 
1918)  for  the  rotating-disk  and  Sagnac  type  of  ex¬ 
periments,  and  the  former  represents  a  plane 
gravitational  wave  which  carries  energy-momen¬ 
tum.  There  are  other  interesting  exact  solutions 
one  of  which  seems  to  describe  cylindrical  gravi¬ 
tational  waves.5  A  general  iteration  solution  (1) 
with  the  boundary  -oi.ditions  that  at  the  observa¬ 
tion  point  g tl„  goes  over  into  a  Lorentz  form  ?)pt, 
is  given  by 

UciV,  Vva^P i>  ]  (iff 

■4  fbh.  x<P-Vic,<P\a- +  (12) 


which,  in  a  case  where  c>0°  =  v  is  the  largest  com¬ 
ponent,  can  be  approximated,  valid  through  sec¬ 
ond  order,  by  the  simpler  form 

AV  =  1„,U  *  2</>  +  2c 2)  -  44%, .  (13) 

Note  that  this  approximation  differs  already  by 
in/,',,  from  Einstein’s  theory.  Other  higher 
order  terms  are  of  the  form  -  *i/00i/>3  for  g00  and 
~  8r?ooV’,  c<,  for  g0„. 

Concerning  the  three  classical  tests  of  general 
relativity  we  may  point  out  that  the  solution  (10) 
agrees  with  the  isotropic  Schwarzschild  line  ele¬ 


ment  in  l’.-st  and  second  order  for  gw  and  in  first 
order  for  gu.  Since  the  three  crucial  tests  arc 
insensitive  to  higher  terms  in  these  components 
the  present  theory  explains  the  three  classical 
effects  just  as  well  as  Einstein’s  theory. 

With  regard  to  the  experiment  on  the  isotropy 
of  inertia  the  present  theory  seems  to  be  in  a  bet¬ 
ter  position  than  Einstein’s  theory  because  our 
line  element  is,  under  (7),  isotropic.  It  may  be 
argued  that  live  principle  of  equivalence  really 
implies  isotropy  of  space  because  the  gravitation- 
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al  energy  is  by  definition  only  a  function  of  posi¬ 
tion  whereas  the  inertial  energy  could  depend  on 
direction  if  the  velocity  of  light  were  different  in 
different  directions.  Thus  the  equality  of  inertial 
and  gravitational  energy  (mass)  seems  to  imply 
the  isotropy  of  the  velocity  of  light,  lienee  the 
■'.utoniatically  isotropic  functional  solution  (10) 
could  be  said  to  represent  more  directly  the  ex¬ 
periment  on  the  isotropy  of  inertia.  For  this 
interpretation  to  be  tenable  one  must,  however, 
assume  that  our  usual  definition  of  inertia  corre¬ 
sponds  to  the  locally  special  relativistic  coordi¬ 
nates  implied  by  (1)  and  (7).  (Although  the  gener¬ 
al  covariance  is  not  harmed  by  (7)  and  (1)  they 
nevertheless  seem  to  imply  that  the  natural  laws 
are  restricted  to  those  exhibiting  locally  Lorontz- 
invariant  field-theory  forms  in  the  coordinates 
described. ] 

For  the  time  delay  of  radar  signals  echoed  from 
Venus  there  seems  to  be  a  curious  difference  in 
the  predictions.  This  difference,  however,  does 
not  arise  from  a  difference  in  the  terms  or  ap¬ 
proximations  of  the  line  elements  but  rather  from 
a  difference  in  the  interpretations  of  measure  in 
the  two  theories,4  In  the  usual  theory  one  identi¬ 
fies  the  zeroth-ordor  special  relativistic  time 
delay  from  the  standard  nonisotropic  solution  of 
Einstein’s  theory*  whereas  under  (7)  there  are 
only  isotropic  solutions  and  our  correspondence 
procedure  seems  to  imply  that  the  zeroth  order 
delay  should  be  taken  from  our  isotropic  solution 
(10).  It  is  then  found  that  the  time  delay  seems 
to  be  4 M/c3  -20  nsec  larger  than  the  value  given 
by  Ross  and  Schiff.J  We  believe  the  problem  is 
woithy  of  further  clarification,  especially  since 
it  is  not  clear  whether  the  actual  experiment  uses 
the  nonisotropic  line  clement  for  correspondence 
or  implicitly  the  isotropic  one.  The  experiment 
at  present  seems  to  give  -5  —  10  nsec  above  the 
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Ross-Schiff  value  but  probable  errors  are  of  the 
order  of  20  nsec.  The  situation  should,  how  eve r , 
be  improved  greatly  within  few  years  to  an  errem 
of  3  nsec. 

From  a  study  of  the  vanishing  covariant  diver¬ 
gence  of  (5)  it  is  concluded  that,  in  this  theory, 
in  order  to  obtain  the  geodesic  equations  of  mo¬ 
tion  (8),  one  must  have  in  the  limit  of  static 
fields  the  condition  6 m  +m6</>=0,  which  appears 
to  be  just  the  expression  of  variation  of  rest 
mass  under  the  principle  of  equivalence  (gravita- 
tion  of  energy).  Thus  the  interpretation  given 
hei  e  to  the  spin-2  field  as  the  expression  of  grav¬ 
ity  proper  seems  justifiable.  (In  Einstein’s  theo¬ 
ry  the  corresponding  condition  is7  6 m  =0.)  If  un¬ 
der  a  more  extensive  analysis  this  interpretation 
of  ours  turns  out  to  be  impossible  the  field 
will  have  to  be  interpreted  as  an  extraneous  field 
in  interaction  with  the  gravitational  field  £„,,(*). 

The  present  theory  is  also  derivable  from  a 
general  action  principle. 
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APPENDIX  2 


PLANE  GRAVITATIONAL  WAVE.' 


Let  600  g33  =  Sj!  =  -  1  +  4Q,  g22  1  -  4Q,  g12  =  4R 

where  R  and  Q  are  functions  of  t  -  z  alone.  This  is  the  usual  plane- 
wave  solution  in  the  linearized  Einstein  theory.  It  is  also  a  first 
order  solution  for  the  new  theory.  Let  us  now  investigate  what  happens 
in  second  order  and  beyond  in  Einstein's  theory.  In  second  order  we  let 


=  -  1  +  4XQ  +  X2U 

(1) 

=  -  1  -  4XQ  +  X2V 

(2) 

=  4XR  +  X2W 

(3) 

where  U,  V,  W  are  quadratic  forms  in  R  and  Q  and  X  is  an  expansion 
parameter  later  to  be  set  equal  to  unity.  The  Einsteinian  equations 

=  0  (4) 

plus  the  wave  conditions  lead  to  the  single  equation 


(U  +  V)"  +  32(QQ"  +  RR")  +  16  (Q ' 2  +  R'2)  =  0  (5) 


where  prime  denotes  partial  derivative  relative  to  t  or  z.  However 
there  exists  no  solution  to  (5)  in  the  form  U  +  V  =  aQ2  +  gQR  +  yR2 


12  - 


except  for  m  :•  '  *-'al  cases  of  R'  =  Q*  -  0  (space  is  flat)  and  R  = 

iQ  (field  is  v>'C  .:eal).  In  the  latter  case  there  is  the  additional 
objection  that  R  -r,c  '  are  not  linearly  independent.  Other  arguments 
show  that  no  function 


U  +  V  =  F (R,Q)  (6) 

where  R  and  Q  are  linearly  independent  can  satisfy  Equation  (5).  Of 
course,  if  the  requirement  of  two  linearly  independent  functions  is 
dropped  (spin  becomes  zero)  there  exists  the  possibility  -g  = 

x  ) 

“S22  e  ,  U  +  V  =  -  16Q2  -  16K2  •  8K,  but  the  linearized  equations 

already  show  that  K  =  At  +  Bz  +  C.  Then  from  (5)  Q  itself  becomes 
prescribed  (Q  contains  no  information)  . 

In  the  new  theory  this  same  problem  goes  as  follows:  The  metric 
is 


Syv  ~  npv  e 


2(4  -  24) 


(7) 


4  = 


Q  R 

R  -Q 


(8) 


,  u  n  ~ 

where  (^  =  ^  =  0,  =  0. 


Expansion  of  (7)  yields 
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,:i  1  ”  “  cosh  4ft  +  ^  sinh  4ft 

?22  “  -  cosh  4ft  -  ^  sinh  4ft 
R  .  , 

gj2  =  sinh  4ft 
g00  g33=1 


(9) 

(10) 
(11) 
0  2) 


1/2 

where  ft  *  (R2  +  Q2)  .  This  metric  i;  proven  to  be  an  exact  solution 

the  field  equations  of  the  new  theory.  The  wave  equation  is  indeed 
satisfied,  as  stated,  and  the  stress-energy  tensor  is  the  canonical 
stress-energy  tensor  of  the  Lagrangian 


t  -  2  J  ,  1 

Lf  ‘  *°  3X*a  +  2  3  * 


(13) 
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GENERAL  STATIC  CURVED  SOLUTIONS  OF  EINSTEIN'S  FIXATIONS  FOR 
SPHERICALLY  SYMMETRIC  METRICS* 
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For  a  general,  static,  spherically  syn- 

inetric  co-ordinate  system  we  find  the  general 

curved  solutions  of  the  Einstein  equations 
.  a  n  .  .  n 

Op  —  u  cuiu  -•  Avp .  r  j  uni  cucse  gcnei ax  su- 
lut  j  ms  the  well  known  standard  solutions  are 
easily  generated. 
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INTRODUCTION 


In  a  c' .  soherically  symmetric  coordinate  system  there  are 
a  number  o'.  *.  "  ’  known  particular  solutions  of  Einstein's  equations 
corresponding  ‘ c'.'^eific  choices  of  the  energy-momentum  tensor. 
However ,  the  problem  of  finding  the  general  solution  of  the  Einstein 
equations  corresponding  to  a  particular  matter  distribution  does 
not  seem  to  have  been  considered.  It  could  be  argued  that  in  view 
of  Birklioff 's  theorem*  this  problem  would  he  of  marginal  interest. 
This  is  not  the  case,  however.  On  the  one  hand  the  various  proofs 
of  the  Birklioff  theorem  do  not  specify  a  method  for  constructing 
the  most  general  form  of  the  metric.  Secordly,  it  is  of  mathematical 
interest  to  find  general  solutions  for  any  system  of  non-linear  dif¬ 
ferential  equations  as  the  existence  proof  might  iriss  some,  parti¬ 
cular  solutions.  Finally,  although  the  general  solution  may  itself 

uc  ci  i_  i  o  1 1  &  j  uj  :■> uine  pdi  Cuj  di  bOluLiOii  ,  Ll i u i  v  j. i iu  cl  iu  i ui  i 

guarantee  that  all  solutions,  for  a  given  matter  distribution  will, 
give  the  identical  predictions  of  physical  phenomena.  This  final 
point  is  discussed  elsewhere^  where  we  show  that  certain  transforms 
ol  the  Schwarzchild  solution  do  not  seem  to  lead,  unarubigously ,  to 
the  same  physical  predictions. 


Below,  we  solve  the  Einstein  equations  for  a  spherically  sym¬ 
metric  coordinate  system.  We  consider  both  the  vacuum  equations 
and  the  equations  with  the  cosmical  constant.  From  these  general 
solutions  certain  well  known  solutions  are  easily  generated. 

All  mathematical  statements  in  this  paper  have  been  verified 
and  intermediate  steps  simplified  using  the  symbolic  manipulation 
system,  MACSYMA,  developed  by  the  Math  Lab  group  at  Tro-ject  MAC, 
M.I.T. 
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1  •  VACUUM  EQUAT IONS 

Consider  the  spherical  coordinate  system  in  which  the  V 
element  may  be  written  in  the  form14 


ds?  =  -  A  dr2  -  F (dO2  +  sin26d^2)  +  D  dt2  ,  (1.1) 

where  A,  F,  D  are  radial  functions  of  class  C2 .  We  wish  to  find 
the  most  general  relation  between  A,  F  and  D  so  that  (1.1)  satis¬ 
fies  the  Einstein  vacuum  field  equations.  We  shall  not  be  inte¬ 
rested  in  cases  in  which  A,  F  or  D  themselves  vanish  everywhere. 

It  can  be  shown  that  for  the  metric  (1.1)  the  only  non  vanish¬ 
ing  components  of  the  Einstein  tensor  are 


G3  =  D MAF  -  (?')■  L -  2FF'))' 
4AF2)j 


(1.2) 


G2  =  G3  =  jl2J'-(2FF/ '  ~  (E/)?)-A/FF/]+])[A(])/FF ,H-2D"F2)-A,I)/F  l-AF(p')2 


2  3 


/♦A2!'2])2 


(1.3) 


G1'  =  4A2!’2  +  ?A'FF'  +  Af(F')2  -  /,FF"1 
“  4A2F2 


(1.4) 


These  components  set  to  zero  yield  the  differential  equati 


quations  we 


wish  to  solve. 
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We  approach  this  system  of  equations  by 
may  be  solved  for  A  as 


.ting  that  Gj  = 


F '  (DF '  +  2D‘  F) 
4DF 


(1.5) 


We  substitute  this  into  (1.3)  and  (3. A)  (set  to  zero)  to  find 


„2  „3  D'F(2D'FF"  -  3T)'  (F7) 2  -  2D"FF') 

^  _  _  o 


'2  3 


(F,)2(DF/  -I-  2D'F) 


2 


(1.6) 


and 


g4  2D(2D,FP"  -  3D'(F/)?  -  2D'  'FF')  = 
**  (F' )  (DF '  +  2D'F)2 


(1  7) 


It  is  clear  that  if  D  /  0  and  F  /  0  the  most  general  way  we  may 
satisfy  (1.6)  and  (1.7)  simultaneously  is  by  requiring 


2D'FF"  -  3D'(F')2  -  2D"FF'  =  0 


This  differential  equation  is  easily  integrated  to  yield 


D  =  Kj  -  K21 


.-3/2 


where  Kj  and  K2  are  constants. 


-h.-  - 


ll  Mild 


(1.8) 


(1.9) 


. — ^ - .......  .. 


In  (1.6)  and  (1.7)  we  note  that  F'  ^  0  and  .  "D'F  i  0.  These 

restrictions  on  the  functions  arc,  from  ;,o..b  Lng  more  than 

the  requirement  A  4-  0  which  we  stipulated  above. 


We  now  substitute  (1.9)  into  (1.5)  to  obtain 


'  \  2 


A  “  — 


K,(F') 


'I”  (Kj  -  K2F_1/^ 


(1.10) 


From  (1.9)  and  (1.10)  in  adlition  to  a  trivial  coordinate  transforma¬ 
tion  we  then  have-* 


TUI- OREM:  The  metric  tensor  defined  bv  the  line  element 


ds2 - - - t-tt-  dr2  -  F (dO2  1  sin20d£2) 

4F(.l  -  KF~1/?) 


(F')? 


+  (1  -  KF~1/?)dt2  ,  (1.11) 


saj_isfie£  tin-  Einstein  vacuum  eqnati  ons  identically  where  F  is  an 
arbi  trary  function  of  r_  of  class  C2. 


We  may  transform  (] 
F  ->  B2r2  where  B  -  B(r) 


.il)  into  a  more  intuitive  form  by  letting 
is  arbitrary.  With  this  new  definition  it 


is  easily  seen  that 


„  _  .1  OLcIll! 


i  -  -K- 
Br 


dr2  -  B?(r?dO?  +  r?sin20dsS?)  b (1  -  — )dt2 

Br- 


(1.12) 
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The  non-isotr. ;  4  •  ^  -bwarz child  solution6  may  now  be  obtained  by  set 
ting  B  =  1  to  £.■ 


ds2  =  ~ 


1  - 


V-  dr2  -  r2 (dO2  +  sin20d({)2)  +  (1  -  -)dt2 


(1.13) 


From  U-ll;  we  may  also  obtain  Fock's  solution7.  By  choosing 
F  =  (r  +  K,  _;2  we  find 


,  K 
r  +  — 


K 

ds2  =  ~  dr 2  -  (r  +  |)2(d02  +  sin20d<;,2)  +  §  dt2.  (1.1/,) 


r  "2 


r  +  f 


From  (1.11)  or  (1.12)  we  may  find  the  isotropic  spherically 
symmetric  iine  element  dv  settinp 


(F,)2 


4F(1  -  KF"1^2)  r2 


(1.15) 


This  differential  equation  may  easily  be  solved  for  F  and  we  find 


F  = 


<*1  + 

«? 


(1.16) 


whore  Kj  is  an  integration  constant.  Thus,  the  isotropic  line  element 
takes  the  form 
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k  ,*• 


ds2  =  -  (1  +  )  (dr2  +  r2d02  +  r2sin26d<{>2)  + 


dt2  (1.17) 


in  accordance  with  the  well  known  result.8 


2.  EQUATIONS  WITH  THE  COSMTCAL  TERM 

For  the  metric  (1.1)  it  follows  from  (1.2),  (1.3)  and  (1.4) 
that  the  component  form  of  the  Einstein  equation 


G1  =  A6* 
J  j 


(2.1) 


where  A  is  ti.e  "cosmical  constant"  is 


r1  =  D(4AF  -  (F')'’)  -  2FF/I)/  _ 
i  -  A 

4AF2D 


(2.2) 


n2  .  r3  .  D^A(-lF,/-(F/)^)-A/rF,HDfA(D,rF/-f2D//F2)-A/D/FZl-AV2(nM2  . 

23  ,990  -  -  A  (2.3) 

4A2F2D2 


4A2F2  +  2A/FF/  +  A((F/)2  -  4FF' ' )  _ 
4A2F2 


(2.4) 


As  we  did  fo  the  more  simple  case  we  solve  (2.2)  for  A  to  find 


F' 

4FD 


])];/  •;  2D'F 


) 


A  ~  - 


FA  -  1 


(2.5) 
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We  substitute  this  into  (2.3)  end  (2. A).  In  both  cases  we  then 
subtract  (the  right  hand  side)  from  the  left  to  obtain  differen¬ 
tial  equations  which  turn  out  to  be  similar  to  each  ether.  This 
procedure  gives 


2.5 


-  \  -  .2p{  (2D/F_2F2j:Il(r/)  (r,)?'-2D/  /F2F/)X-2n/VF/  /-j-3D/  (F')  ^-!-2D/  'FF 

F'(DF'  +  2D'F)? 


(2.6) 

and 


-  A 


-  P  1*1(20  F^I  +D(F/)  -D/F(F/ )  ^-2D/ /F2F ' )  X— 21) ' FF '  ,-L'm '  /f ' \  2 


,-h3D/(F/)  42D'1 


(F')  (DF'  +  2D'F)2 


(2.7) 


It  is  clear  that  to  solve  (2.6)  and  (2.7;  we  need  to  find  the 
relation  between  F  and  D  ;  iiich  causes  the  expression  in  curly  brae 

kets  to  vanish.  For  1  =  0  we  know  that  the  relation  we  seek  must 
reduce  to 


D  = 


K, 


hence  we  now  seek  a  function  N (F)  such  that 


»  =  Kj  -  K2F_1/2  +  AN(F) 
satisfies  (2.6)  and  (2.7). 


(2.8) 


(2.9) 


FF'} 


Substitutir  .>  (2.9)  into  (2.6)  results,  after  some  laborious 


calculation,  is  th '  f.  '.lowing  differential  equations  tor  N: 


N  -  N'F  =  0 
3N'  -  K,  =  0 


(2.10) 


These  are  easily  integrated  to  find 


N  = 


Thus,  (2.9)  becomes 


(2.11) 


n  «*  v  n 

l ' 


v  V-'!2 
1 


(2.12) 


With  this  expression  for  D  we  find 

K  (F')2 

A  - - - - - - 

4(Kj  -  K2F"1/2-  |  KjF)F 


(2.13) 


From  (2.13)  and  (2.14)  and  a  trivial  coordinate  transforma¬ 
tion  to  absorb  an  integration  constant  we  have'’ 

THEOREM :  The  metric  defined  by 


ds2  = 


(FT 


4*F(1-KF“1/2-  ~  F) 


dr?  -  F (do2  +  sin20d$2)+(l *-  KF-1/?-  ~F)dt2 


(2.14) 
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satisfies 


Gl  X{lJ 

v  v 


(2.15) 


for  any  function  F(r)  of  class  C2, 


Lot  us  now  choose  F  =  r2  in  (2.1/).  With  this  choice  we 


find 


ds2  =  - 


dr2 


1  .  K  .  X  ~  -  r2(<l“2  sin?0d(f 2)  +  (1  -  f  -A  r2)dt2.  (2,16) 

r  3 


Tins  is  the  Kottler  solution9  and  represents  one  of  the  very  few 
known  (if  indeed  other,  have  been  found)  solutions  of  (2.1). 

We  now  derive  an  expression  for  the  isotropic  solution  cor¬ 
responding  to  (2.14)  by  setting 


(F')' 


4F(1  -  KF“l/?-  A  F) 


(2.17) 


We  see  that  F  must  satisfy 


/, 


dF 


F(1  KF  1/2-  ^  F)1/2 
3  ’ 


-  =  In 


c*r 


(2.18) 
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"here  C  U  an  Integration  constant.  This  does  not  appear 
^  he  expressible  in  dosed  ion,  However,  one  ear.  obtain  expres¬ 
sion  for  particular  cases  of  interest.  By  setting  K  -  0  in  (2.18) 
we  may  perform  the  integration  and  solve  for  F  to  find 


Then 


12  CR2 
(XR2  +  C)2 


(2.19) 


ds2  =  - _ 12  C 


(AK2  +  C)2  Wr2  +  r2<i02  +  r^''20d^)  +  (1 _ ACr£ 


(AR2  +  C) 


-)dt2 


Here  we  have  the  isotropic  st-nt-iV  f 

[UC  static  form  of  the  deSirt^*-  .  • 

aeMtter  solution.  Othe 

00  uLlons  of  possible  cosmological  interest  could  be  found  by 
Performing  a  series  expansion  on  the  left  Hand  side  of  „.le). 

«e  Have  not  considered  such  cases  at  the  present  time. 


(2.20) 


M 
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APPENDIX  4 

UNIFORM  ACCELERATION 


Consider  the  Cartesian-like  coordinate  system  defined  by 


ds2  =  -  F(dx?  +  dy2)  -  Hdz2  +  K  dt2 


(1.1) 


where  F,  H  and  K  are  function  of  z.  He  wish  to  find  the  general 
relations  between  these  functionals  so  that  (1.1)  satisfies  Einstein's 
Vacuum  Equation: 

For  (1.1)  we  find  the  components  of  the  Einstein  tensor  (set  to 
zero)  to  be 


o1  -  r2  _ 


!  I'  ^11  f  1P?1>  '  t  i 


....  .  PPti.li  rPi.il.  I\ 

1  l  4V  '  14  J  A  *1  IS.  J 


ft*  ?  \  ? 


+  K2  (U (2FF '  '  -  (F*)2)  -  FF'H ' )  }  /  4F2112K2  =  0  (1.2) 


. 3  _  _  F'(2FK'  +  KF') 
3  4F2HK 


=  0 


(1.3) 


G4  = 
4  . 


H  (F[2F2F '  '+F(F'  )21  +  F2(2IT'  (F ' )?)  -  F2  (F'  )  2  )-2F3F '  U ' 

4F4il2 


We  may  satisfy  (1.3)  by  requiring 


(1.4) 


2FK'  +  KF'  =  0 


(1.5) 


This  diff aquation  may  be  integrated  and  we  find 


F  =  Cj/K2  (1.6) 

where  Cj  is  a  constant.  This  relation  implies  that  (1.2)  and 
(1.4)  become 


G 


1 

1 


11(2KK"  -  5(K')2)  -  KH'K' 
4H2K2 


(1.7) 


G 


*4 

4 


1I(2KK//  -  5(K7)2)  -  KH'K' 
H2K2 


(1.8) 


Hence,  by  solving  (1.7)  or  (1.8)  we  will  have  a  general  solution 
of  the  Einstein  Vacuum  Equation.  These  are  integrable  and  we  find 

CIIUL 

GO2 

H  =  C,  -  (1.9) 

2  K5 

satisfies  (1.8)  and  (1.8).  Thus,  the  metric  tensor  components 
defined  by 


ds2 


—  (Jx2  +  dy2)  -  C.  -----  dz2  +  K  dt2 
K2  2  K5 


(1.10) 


is  the  general  static  solution  of  G^  =  0  for  any  function  K(z)  4  0. 


—  - . 
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On  the  other  hand,  we  may  satisfy  (1.3)  by  choosing  F'  =  0. 
With  this  choice  it  is  found  that 


=  _  2KH  K"  +  -K  H/K/  - 
4i(2K 


=  0 


(1.11) 


and 


«*  0 


The  solution  of  (1.11)  is 


(1.12) 


Cj(K')2 


We  thus  find 


ds2  -  -  (dx2  +  dy2)  -  C,  dt2  +  K  dt2 


also  satisfies  G*  =  0 

j 


(1.13) 


(1.14) 


llouevcr,  it  is  not  difficult  a  prove  that  for  the  metric  (1.14), 
all  components  of  the  Rieman„  Christoff cl  Tensor  vanish  identically. 
Hence,  the  metric  (1.14)  is  flat.  It  can  also  he  shown  that  (1.14)  is 
the  only  flat  metric  belonginB  to  the  general  class  (1.1). 

From  (1.10)  we  m„y  seek  the  general  isotropic  line  element  by 
setting 


c,,(kV- 

K2 


(1.15) 


Then, 


K'  =  C  K372  (1.16) 

or 

K  13  (Cx  +  D)2  (1.17) 

Hence 

ds2  -  -  (Cx  +  D)  4(dx2  +  dy2  +  dz2  ) 

+  (Cx  +  P)?  dt2  (1.18) 


is  the  isotropic  line  element. 
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APPENDIX  5 

* 

NEW  SOLUTIONS  F  0  R  E  I  N  S  T  E  i  :  T  11  E  0  R  Y 

11  us  on  in  Yilnuiz 

Perception  Technology  Corporation ,  Winchester ,  Nassacl'usctts  01890 

(Analytical  Subject  Index  01.0) 

ABSTRACT 

A  general  solution  for  Einstein' s  theory  of  gravitation  corres¬ 
ponding  to  spherical  symmetry  is  presented .  It  is  shown  that  the 
usual  Schwarchi Id  metrics  are  members  of  this  general  solution. 

Some  of  the  new  solutions  implied  seem  to  make  little  or  no  physical 
sense.  One  that  at  the  outset  seems  to  be  physically  unobjectionable 
and  interesting  for  possible  astrophysical  applications  is  described 
in  some  detail. 

By  the  use  of  a  well-known  method  of  tensor  transformation  it  is 
proven  that  i.hc  line  clement 


ds 2  -  fc'-Qt^  - 


,2m 


(d  (r  +  sin/ Odif' )  -  f 


1  2m  ,? 
[l-fJ 


dr * 


(1) 


where  f  is  a  twice  differentiable  function  of  r  is  a  general  solution 
of  Einstein's  field  equations  R  =  0.  rite  usual  Schwarzchild  solu- 
tionsi  art?  members  of  this  general  solution  for ^ 


TO  BE  SUBMITTED  FOR  PUBLICATION 


32  - 


f  =  1  -  2<j>  ,  £  =  1,  -  2<$(1  -t  p/2; 


(2) 


where  4>  =  m/r  so  Ih at  for  large  r  the  first  orae-  Etonian  corres¬ 
pondence  is  satisfied.  Similarly  f  =  (1  -  P)/(l  +  p)  yields  the 
well-known  Fork  solution 7  satisfying  the  correspondence  condition  in 
first  order.  The  solution  (1)  is  new  and  offers  enormous  flexibility 
as  we  shall  see  }>clow.  It  was  originally  conjectured  during  the  per¬ 
formance  of  an  ARPA  contract  directed  to  study  the  solutions  of  au¬ 
thor's  new  theory  of  gravitation-1 ,  and  later  explicitly  verified  by 
the  use  of  the  MACSYMA  system  of  Mathematic' s  Laboratory ,  Project 
MAC,  M.l.T.  MACSYMA  was  also  of  definite  value  in  the  formulation  of 
the  conjecture  itself. 

Actually  it  is  somewhat  disturbing  to  find  that  f  can  be  so 
general  as  in  (1) ,  because  we  seem  to  be  able  to  set 


f  -  >;  a  pn  (3) 

n 

and  produce  an  arbitrary  number  of  Schwarychild  type  singularities , 
or  choose, 


(4) 


and  obtain  a  solution  with  no  singularity  at  all,  except  at  r  -  0. 

Note  that  some  of  the  possibilities  contained  in  (3)  might  be 
problematic  for  the  experimental  aspects  of  the  theory.  For  example, 
Inj  setting 
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element  for  all  experimental  purposes,  hcr,c'~  !t  reproduces  all  four 
of  the  experimental  tests,  b)  In  the  limit  r '  :,low  motion  of  the 
test  body  the  geodesic  equations  of  motion  redan  *  r 

d?x 

- £  =  {  °)  =  -  3  {  (9) 

dt 2  y°  °  y 

which  arc  exactly  the  Newtonian  equations .  c)  In  the  limit  of  slow 
motion  of  the  source  the  field  equations  reduce  to 

V2<t>  =  -  4va  (10) 

which  is  exactly  the  Poisson  aquation.  By  contrast  other  line  ele¬ 
ments  reduce,  to  these  Newtonian  limits  in  an  approximate  sense  only, 
d)  It  is  the  only  solution  of  Einstein's  equations  which  apply  both 
*-  ~  4-u  r ..  4. «  —  ,*  ■*.  -  — i  j  >  -■  —  r-  r-  -■ —  zzz’cy.s  of  zl  ziz.ee  f  le  t  ±S.f'Z'tizte  *,\i  t 

out  a  change  in  form,  c)  It  is  harmonic  in  the  sense  that  if  the 
metric  is  written  in  the  form,  ds 2  -  Dc2dt2  -  B(r2dQ2  +  r2sin2Qd$2 )  - 
Adr2 ,  that  is  as  a  factor  to  special  relativity ,  then  rd^(\'-g  gVV )  - 

8  (>/DB2A  A  1)  -  0.  f)  It  has  no  singulari  ty  anyv/hcrc  except  at 
v 

r  •>  0  where  i  t  reduces  to 

ds2  «  ~(2m)2(dQ2  ■ h  sinHdv2)  (11) 

The  latter  represents  the  surface  of  an  impenetrable  sphere  of  radius 
2m.  This  spherical  surface  seems  to  be  suspiciously  similar  to  the 
Schwurzchi Id  surface  but  note  the  difference  that  (8)  has  an  essen¬ 


tial  singularity  at  the  surface,  not  an  oidinary  one.  Fu:  thermorc 


there  is  really  no  interior  of  this  sjihere  as  r  starts  out  at  the 
surface.  The  only  conceivable  objection  to  (8)  might  be  that,  like 
all  othei  solutions  of  Einstein's  theory ,  it  leads  to  a  zero  gravita¬ 
tion,  ll  stress-energy  tensor  £  and  thereby  seems  to  violate  the  New- 
torn  an  corrcsp  ndcncc  because  we  arc  missing'* 

tv  "■  ~V’  3  *  +  2  0  *  V  (12) 

As  to  a  possible  application  of  (8)  other  than  the  usual  astro¬ 
nomical  calcuJ ations  one  naturally  thinks  of  the  gravitational  col¬ 
lapse.  It  can  be  seen  that  for  (8)  there  does  not  exist  the  concept 
of  a  black  hole  in  the  usual  sense.  The  collapsed  body  would  here 
remain,  at  all  limes,  as  a  naked  object.  Furthermore  there  is  no 
singularity  at  finite  x  and  the  usual  reversal  of  the  signature  of 
the  metric  dons  not  Lake  place. 

It.  should  be  interesting  to  repeat  the  usual  calculations  of 
gravi  tat.i  onal  collapse  and  the  dynamics  of  black  holes  using  some  of 
the  new  members  o  the  general  solution  (1),  including  the  unusual 
member  (8) . 

It  is  a  plc.isuuc  to  thank  Dr.  Richard  Pavel lc  for  tcclmical 
help  and  discussions . 

Inis  ri 'search  was  supported  by  the  Advanced  Research  Projects 
Agency  oi  the  Department  of  Defense  vndci  Contract  No.  DAHC15  73  C  0369. 
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That  the  conjectural  line  element  (1)  is  a  general  curved -s; pace  solu¬ 
tion  Oi  hinstein' s  equations  was  proven  bj  <.  r  '•  er  transformation  on 
the  metric.  There  are  also  purely  flat-spa _  .  solutions  of  the  form 
ds?  ~  c?dt  -  F2  (d 02  +  sin20d$2)  -  F,2dr2  which  c~c  ot  discussed  in 
this  paper.  These  solutions  arc  checked  by  KAiCSYKA,  a  symbolic  mani¬ 
pulation  system,  at  Project  MAC,  MIT.  The  author  is  grateful  to 
Dr.  Joel  Moses  and  Mr.  David  Grahcl  of  Project  MAC  for  making  this 
verification  ,x>ssible.  A  direct  proof  of  (1)  is  found  by  Dr.  diehard 
Pavel lc  Oi  Perception  Technology  Corporation  by  using  the  extensive 
manipulative  pciwcr  of  MACSYMA  system.  Whether  these  are  the  only  pos¬ 
sible  gcnei  a  1  solut  ions  of  the  theory  in  spherical  coordinates  has  not 
yet  been  established,  although  a  reasoning  leased  on  Dirkhoffs  theorem 
seems  to  indicate  that  there  are  no  others. 

1.  Any  text  hook  on  General  Relativity. 

2.  V.  Fock,  "The  Theory  of  Space  Time  and  Gravitation" ,  p.  215, 
Pergamon  Press,  (1964). 

3.  II.  Yilmaz,  Phys.  Rev.  Lettrs.  20,  1399  (1971);  II  Nuovo  Cimento, 
10D,  79  (1972) 

<7.  hoi  a  method  of  calculation  see  for  example  S.  Weinberg,  ‘Gravi¬ 
tation  and  Cosmology" ,  Wiley  (i  Sons  Publishing  Co  (1972). 

5.  If.  Yilmaz,  Lcttcrc,  Al  Nuovo  Cimento,  7 ,  337  (1973). 
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APPENDIX  6 

INDTC.'.AL  MANIPULATION  IN  THE  HACSYMA  SYSTEM 


The  HACSYMA  system  consists  of  over  250  functions  for  the  manipula¬ 
tion  of  mathematical  expressions.  There  are  not  only  functions  for  per¬ 
forming  mathematical  computations  but  also  those  for  many  necessary  mis¬ 
cellaneous  manipulations  such  as  extracting  parts,  simplifying,  input 
output,  etc.  Each  one  of  these  is  made  up  of  one  o'  as  many  as  several 
hundred  programs  (xn  the  case  of  integration)  written  in  the  LISP  lan¬ 
guage.  The  expressions  dealt  with  are  composed  of  operators  and  argu¬ 
ments  to  them  which  are  numbers,  variables,  functions  references,  arrays, 
lists,  and  matrices.  The  mathematical  knowledge  built  for  MACSYMA  is 


wntumflii,  cAp.-.uuiug  as  MOie  aim  muic  piugiam.s  din  wiixLeii  to  pjoviue 

additional  capabilities.  At  present  there  are  over  2000. 


If  one  wishes  to  introduce  a  new  mathematical  entity  into  the  system 
then  he  would  not  only  have  to  cast  it  in  terms  of  the  existing  syntax 
but  also  it  would  be  necessary  to  provide  the  programs  to  perform  the 
desired  manipulations  on  it  as  it  would  be  foreign  to  most  of  the  exist¬ 
ing  programs.  (MACSYMA  does  have  a  facility  for  syntax  extension  but 
currently  this  applies  only  to  operators.)  Perception  Technology  Corpora¬ 
tion  has  embarked  on  an  effort  to  create  some  facility  for  dealing  with 
tensor  manipulation  as  MACSYMA  has  nothing  that  would  provide  a  conve¬ 
nient,  efficient,  and  workable  mechanism  for  indicia]  types  of  manipula¬ 


tions  . 


way  it  is  written  bit  by 


Since  a  tensor  is  defined  not  by  the 
certain  properties  it  has,  we  represent  them  as  more  general  quanti¬ 
ties  which  we  call  indexed  objects.  They  possess  a  name  and  two  lists 
of  indices,  the  first  list  being  the  subscripts  (covariant  indices) 
and  the  second  list  being  the  superscripts  (contravariant  indices). 
These  indexed  objects  are  represented  as  functions  of  two  arguments, 
which  are  the  above  mentioned  lists.  If  either  set  of  indices  is 
absent  it  is  represented  by  an  empty  list.  Thus  6^  and  A±ik  would  be 
represented  as  6([i],[j])  and  A([  ],[i,j,k]). 


Ordinary  differentiation  of  an  indexed  object  (say  E)  with  respect 
to  the  coordinate  Xk  would  be  obtained  by  the  command  DIFF (E ,k) .  Since 
MACSYMA  wouldn't  know  that  E  depends  on  Xk  it  would  give  0  but  the  dif¬ 
ferentiation  program  has  been  modified  to  assume  that  all  indexed  ob¬ 
jects  depend  on  all  coordinates.  If  an  indexed  object  is  independent 
of  all  coorcinates  this  could  be  stated  by  the  command  DECLARE 
(E, CONSTANT)  and  subsequent  differentiation  of  E  would  yield  0.  The 
derivative  of  an  indexed  object  causes  the  coordinate  with  respect  to 
which  differentiation  is  carried  out  to  be  appended  to  the  function 
which  represents  the  indexed  object  in  MACSYMA.  These  indices  which 
denote  differentiation  with  respect  to  a  coordinate  are  sorted  in  alpha¬ 
betic  order  to  take  advantage  of  the  commutativity  of  differentiation 
in  order  to  simplify  expressions  composed  of  those  indexed  objects. 


At  present  two  features  have  been  implemented  -  symmetries  and 
contractions.  For  the  former  purpose  thr.  o  tv.  ions  are  provided. 

One  for  declaring  symmetries  (e.g.  G([  ]>[j,.’-j)  is  equal  to  G([  ]»[i»j]))> 
One  for  removing  them,  and  one  for  displaying  them  in  case  the  user 
fevgot  which  ones  he  declared.  For  the  purpose  of  contracting  products 
of  indexed  objects  in  expressions,  four  functions  were  provided.  One 
to  declare  that  some  indexed  object  conti acts  with  another  to  form  a 
third,  another  function  to  remove  these  declarations,  another  to  dis¬ 
play  them,  and  also  a  function  to  perform  the  contractions.  This  latter 
function  is  rather  difficult  to  code  as  it  is  necessary  to  make  several 
passes  over  the  expression  to  make  sure  that  all  possible  contractions 
have  been  utilized.  The  function  also  performs  the  substitutions  im¬ 
plied  by  the  Kronecker  delta  function  if  it  occurs  in  an  expression. 

In  addition  a  function  has  been  written  fo  display  expressions  containing 
indexed  objects  in  more  ns  :ural  notation  than  MACSYMA  ordinarily  would 
because  MACSYMA  assumes  that  A( [i , j ] , [k] ,l,m)  is  a  function  of  four 
arguments,  and  would  display  it  as  such,  whereas  we  would  like  to  see 

k 

it  displayed  as  A. .  ,  . 

ij,lm 

Using  this  equipment  described  above  we  were  easily  able  to  perform 
the  calculation  of  the  Ricci  tensor  to  first  order  given  the  expressions 
for  the  covariant  and  contravar"1' ant  forms  of  the  metric  tensor.  The 
result  in  unsimplified  form  had  16  terms  but  these  were  reduced  to  6  hv 
using  the  contraction  program. 


To  carry  this  calculation  to  second  order  will  result  in  thousands 
of  terms  if  the  result  was  not  simplified  along  the  way.  Although  many 
times  more  complicated  than  the  first  order  case,  it  will  be  possible 
to  do  but  shall  require  more  code  than  is  available  at  present  for 
simplifications . 


APPENDIX  7 


X 

c  = 


EXPONENTIAL  MATRICES 


In  the  New  Theory,  for  particular  physical  problems,  functional 


X  -  '  X2 

expansions  of  the  form  e  =  1  +  X  4-  -r-p  +  ...  where  X  is  a  square  matrix 


are  considered.  Here  X2  is  the  product  matrix.  The  MACSYMA  system  is 
ideal  for  generating  useful  closed  form  expansions  of  this  kind  and  we 
have  been  able  to  find  many  interesting  cases.  The  conve  ‘gencc  of  the 
following  forms  may  be  proved  by  induction.  In  the  following  examples 
the  dot  represents  a  zero  element. 


-'l  +  (A?+B2)  (cosp-1) 
n2 


Asi np  +  BC(cc  sp-1) 
P 


Bsinp  AC(cosp-l) 


where  p 


2  =  a  2 


A/  +  +  C' 


/ 


X  - 


-A 

-B 


-C 


Asinp  BC(cosp-l) 


l  +  (A?+C2) (cosp-1)  Csinp  AB(cosp-l) 


Csinp  +  AB(rosp-l) 


(1) 


Bsinp  AC  (cosp-1) 


l  +  (B?+C2) (cosp-1)  j 

r.  2  ' 


-  .a  ,  . 


(2) 


(3) 


\-R  -S  -T 


/  1  +  —  (cosp-i) 
/  p2 


—  (cosp-1) 
P2 


—  (cosp-1) 
p2 


sm  p 


RS  . 

—  (cosp-1) 
P2 


1  +  —  (cosp-1) 
p2 


—  (cosp-1) 
P2 


sin  p 


--  (cosp-1) 
P2 


(cosp-1) 


1  +  —  (cosp-1) 
P2 


sin  p 


R  . 

\  "  p  Sln  p 


sm  p 


sin  p 


cos  p 


where  p2  =  R2  +  S2  +  T2. 


A  =  \R  Qi 


x  cosh (R)  sinh(R) ' 


^sinh(R;  cosh  (R)  j 


Q 

•  n  ^ 
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X 

e  = 


i > 

e 


C 

e 


D 


Q  R 


l R  -Q, 


t 

I 


cosh(R'  + 


q2)  1  /2  +  Q  sinh(R?  +  Q2) 
(R2  +  Q2)1/2 


f 


R  sinh(R2  +  q' )  1  ^  2 
(R2  -»•  Q2) 


R  slnh(R2  +  Q2) 1 ^ 2 
(R2  +  Q2)1^2 


\ 


cosh(R2  +  Q2)1/2-  ^_Sinh(R2  +  Q2)1/2 

(R2  +  Q2)  ' 


(8) 


